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Perturbative gadgets were originally introduced to generate effective fc-local interactions 
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1 Introduction 

The anyon based quantum computer, originally proposed by Kitaev [10U161 IS], is a promising 
approach to the realization of fault-tolerant quantum computing. Its main feature is the fact 
that quantum information is encoded in non-local observables related to topological invariants; 
such information is robust against local errors. The physical substrate of these schemes is a 
system with anyonic excitations above a degenerate ground space. Ideally, all computations 
are performed by successive creation, braiding and fusion of anyons, resulting in non-trivial 
actions on the ground space. 

A number of systems exhibiting the required topological ground state degeneracy and 
anyonic excitations are known: These include Kitaev's well-known Z2-toric code and its gen- 
eralization to arbitrary finite groups [TU], Kitaev's model on the honeycomb lattice [2] and 
Levin and Wen's string-net models 12 . The merits of a given proposal may be judged ac- 
cording to the following criteria: (i) the difficulty of realizing the Hamiltonian in a physical 
system and (ii) the computational power of its anyonic excitations. In terms of criterion (i), 
the model [9] is perhaps most exciting, since the terms in the Hamiltonian are 2-local in- 
teractions between neighboring qubits; indeed, corresponding experimental proposals already 
exist [5l 113). Unfortunately, the anyons in this model are not computationally universal. On 
the other hand, Levin and Wen's construction realizes all "doubled" topological phases, giving 
in particular anyons corresponding to all discrete gauge theories and doubled Chcrn-Simons 
theories, but this model generally requires 12-local interactions between qudits and therefore 
seems hard to realize. 
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Perturbation theory is traditionally used to obtain an effective low-energy description 
of a given physical system. Perturbative gadgets turn this procedure around: they were 
introduced in [8] to realize certain "target" Hamiltonians as the low-energy limit of 2-local 
Hamiltonians. This powerful idea has a number of applications [TT1 [JJ [2]; in [17], it was 
observed that perturbative gadgets can generate (arbitrary) fc-local effective Hamiltonians 
from 2-local Hamiltonians. A particularly elegant construction achieving this was recently 
given by Jordan and Farhi [7] based on Bloch's perturbation series [3]. While appealing 
in their generality and simplicity, these results may be of limited use for the realization of 
concrete Hamiltonians because of their large overhead: The number of required qubits and 
interactions scales polynomially in the number of non-trivial summands and factors when the 
target Hamiltonian is expressed in terms of Pauli operators. 

It is natural to assume that gadgets adapted to the internal structure of a given target 
Hamiltonian may lead to more efficient and natural realizations; this is indeed the case for 
certain Hamiltonians that give rise to topological order. Here we focus on Kitaev's model [TU] 
based on a finite group G. Depending on the choice of G, this has the following properties: if 
G is any non-solvable group or a certain semidirect product of cyclic groups of prime order 
(which includes G = S3), then the resulting anyons can be used to realize a universal gate 
set, as shown in a sequence of works [TBI [16j 031 US]- The degrees of freedom in this model 
are (d = |G|)-qudits placed on the edges of a lattice C. The Hamiltonian has the form 

v p 

where for every vertex v, the vertex-term A(y) acts on the incident edges of v, and the 
plaquette-term B{p) acts on the boundary of plaquette p (We will define these operators in 
detail below). In particular, the terms in the Hamiltonian are at most I = max{c?£, dcj-local 
(more explicitly: they act on (C' G ')®'). Here dc and dc are the maximal degrees of vertices in 
the primal and dual lattice, respectively. Concretely, this gives interactions acting on (Cl G l)® 4 
on a rectangular lattice, and interactions acting on (C' G ')® 6 on a hexagonal lattice. 

2 Results and Techniques 

We summarize our results, and sketch the involved constructions. More details will be given 
below; here we merely give a high-level overview of our methods. 

2.1 Quantum double Hamiltonians as effective Hamiltonians 

We construct 

(a) a Hamiltonian H vm t cx which has the sum — ^2 A(v) of the vertex-terms of Hqd as its 
low-energy limit, and consists of terms acting on (C' G ')® 2 £g>C d£ . In a similar way, we can 
obtain a Hamiltonian -ff p i a qucttc which generates — B(j>) and is made of terms acting 

on (d G l)® 2 (g)C^'. 

(b) a Hamiltonian Hf u \\ whose effective low-energy Hamiltonian completely reproduces Hqf>. 
This construction consists of terms acting on (C' G ')® 2 (g> (C')® 3 . 



R. Konig 3 



To clarify the simple structure of these Hamiltonians, we schematically depict the operators 
in iJ vert ex and -ff p i aq uctte which generate a given vertex-/plaquette-operator: 

- + A A + W 

The operators on the right act only on a single edge of the lattice and an auxiliary system 
located at the vertex/center of the plaquette, respectively. This auxiliary system is represented 
by a black circle; for vertices, it is C' G ' <S> C 3 , and C' G ' <S> C 6 for plaquettes. 

The operators in H{ u \\ are obtained by a similar substitution rule, with the difference that 
all auxiliary systems are of the same form C' G ' ® C 6 , and the operators on the right also act 
on auxiliary systems associated with neighboring plaquettes and vertices, respectively: 




- 0+0+0+0OO 

Gray circles indicate that the action on these systems is particularly simple. (In fact, it is 
given by a projection onto a fixed vector |0) £ C 6 and leaves the first part C' G invariant.) 

Comparing the support of the operators in the original Hamiltonian Hqd with those of 
the constructed Hamiltonians, we get the following table: 



lattice 


H QD (original) 


-^plaquette? -^vertex 


Hf u \\ 


rectangular 


ft£|G'|\(g>4 


( C |G|)®2 0C 4 


(Cl«l)® a (g> 


(C 4 )® 3 


honeycomb 


(Cl G 'l)® 6 


(Cl G l)® 2 ®C 6 


(Cl G 'l)® 2 ® 


(C 6 )® 3 



While the dimension of the support is decreased only for \G\ > 8 on the rectangular lattice 
(and \G\ > 3 for the honeycomb lattice), we stress that the terms in our Hamiltonians are 
relatively simple: they involve at most 2 systems related to the group G and are otherwise 
made of hopping-terms. The terms related to the group G are different versions of controlled- 
group multiplication, as explained below. 

Obtaining the vertex-terms. 

Let us sketch the construction of these Hamiltonians for the case of a square lattice. Consider a 
single vertex term A(v); this operator depends on the orientation of the edges of the lattice C. 
For simplicity, assume that v has edges {eg, ■ • ■ , e<j} with arrows pointing away from v. In this 
case, the vertex operator is A(v) = t^t J2 g Ili=o -^-( e i )> w here the operator L 9 _(e^) stands 
for right-multiplication of the degree of freedom on edge e\ by g _1 for g E G. Our aim is to 
find a Hamiltonian H v whose effective low-energy Hamiltonian is equal to A(v). 

We associate with the vertex v an auxiliary system %rv ® %jv = C' G ' ® C 4 with standard 
basis vectors \g)\i), where g e G and i = 0, ...,3. On Hrv, we introduce the projection 
operators 

(T»)«. = \g){g\ R « , geG. 
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We then define the operators 

M*=£z*( e »)®(T^ (3) 

9 

acting on the degree of freedom associated with the edge e v and the auxiliary system T-Lr^ . 
Observe that this is a unitary. It corresponds to controlled-right multiplication of the degree 
of freedom associated with the edge e v by 5 _1 , where g is the value of the auxiliary register. 
Our Hamiltonian is of the form H v = H% + XV V with 

HS = -|* W ><*V ® |0><0|j. , 

V v ='£((MVy®\i + l){i\ I .+h.c,) , (4) 

i=0 

and \^)rv = } J2 a eG Iff) ^ %R°' (Here, h.c. denotes the Hermitian conjugate, and we 
V l G 

omit identities when operators act trivially on subsystems.) We show below that up to a 
constant energy shift, this Hamiltonian reproduces the term —A(v) in 4th order perturbation 
theory in A. 

To obtain the target Hamiltonian — A(v), we apply the described procedure to every 
vertex v. That is, we introduce an auxiliary system T-L v at every vertex, and define Hq = 
J2 V Hq and V = J2 V V V as the sum over vertices of the operators defined by (j4]). This leads 
to claim (jg). 

Obtaining the plaquette-terms 

A similar statement holds for the plaquette-terms B(j>). For concreteness, consider a pla- 
quette whose boundary edges {e^, . . . , e^} are oriented in such a way that they form a 
counterclockwise loop around the plaquette. The plaquette-operator is equal to B(p) = 
2 g3 ...oo=i Ili=o -^-( e ?)' where the the sum is over all 4-tuples (go, . . . ,33) whose product 
53525150 is equal to the identity element 1 G G, and T^_(e p ) = \g~ 1 )(g~ 1 \ e p are projection 
operators acting on the degree of freedom on edge e\. To generate this term in perturbation 
theory, we introduce auxiliary systems Hrp ® Hip = C' G (g C 4 of the same form as before, 
and define the left-multiplication operator 

Z 

on Hrp . We then set 

Mf = X)r_(ef)®(^)«p (5) 
9 

and \^> p )rp = |1) G The plaquette-term —B(p) is then again reproduced in 4th order 

perturbation theory in A from H p = Hq + XV P , with 

flJ = -|* p )<* I, |flp(»|o)<o|j, 

v=x;((jwf) t ®i»+i)<»v+fc. c .) . 

i=0 

The generalization to the target Hamiltonian — -B(p) is straightforward. 



R. Konig 5 



Obtaining the full Hamiltonian Hqjj . 

The previously described naive procedure for generating several vertex- or plaquette-terms 
fails when we simultaneously consider vertex- and plaquette-terms. This is because the opera- 
tors (cf. ((3])) associated with a vertex v generally do not commute with the operators Mj 
(cf. associated with a plaquette p which has v on its boundary. To overcome this problem, 
we introduce additional interactions: For any vertex v, we modify the perturbation (j4]) to 

3 

V v = Y i (W)* ® K + 1><*|J- ® |0><0| jP _(e ? , (8 |0)(0| JP+K) + /i.e.) , 

i=0 

where p_(e") and p+(e£) are the plaquettes separated by the edge e\ (cf. Figure Q] below) ; 
the expressions V p corresponding to plaquette-terms are changed in a similar way. Our per- 
turbative analysis then shows that the resulting perturbation V — J2 V V V + Yl p V P generates 
the desired Hamiltonian Hqjj in the low-energy limit. 

2.2 Perturbative techniques for loop-Hamiltonians 

Our results rely on the following two facts. First, the vertex- and plaquette-terms of the 
quantum double Hamiltonian can be written as products of operators followed by a projection. 
Second, an ordered product of k operators can be obtained in fc-th order perturbation theory 
from a certain hopping-term Hamiltonian. 
The first fact boils down to the identities 

A(v) = {$ V \M£---M£\V") 
B(p) = (* p |M f, ---M 3 p |* p ) 

which can be verified by straightforward computation from the expressions given above. How- 
ever, more insight is gained from a derivation of this statement that uses some fundamental 
objects related to the quantum double Hamiltonian: We first rewrite the vertex-term A(v) as 
a ribbon-operator [10] associated with a closed ribbon (also called Wilson loop) going around 
the vertex v. Operators corresponding to the concatenation of two ribbons may be written as 
a linear combination of products of ribbon-operators corresponding to their parts; this is de- 
scribed by the comultiplication map in the algebra of ribbon-operators. The operators (T? )rv 
and {L 9 + )rp used in the definition of M", are related to a representation of the dual Hopf 
algebra; these provide us with the desired expressions for A(v) and B{p) as products. The 
procedure of writing the Hamiltonian in terms of operators corresponding to closed loops 
(ribbons), and subsequent decomposition of the loops into small segments could extend to 
similar quantum loop models. 

The second fact may be of independent interest. We refer to it as the clock-gadget. 

The Clock-gadget 

This gadget provides a general tool for generating ordered products of operators. It is moti- 
vated by similar hopping-term constructions previously used to prove the QMA-completeness 
of the local Hamiltonian problem [TTJ. Roughly, we have the following statement. 



6 Simplifying quantum double Hamiltonians using perturbative gadgets 



Theorem 1. (Informal version) Let Hq be a Hamiltonian on a bipartite system %$ ® Hi 
of the form Hq = —Pq, where Pq = Tq (g> |0)(0|j is the projection onto the ground space and 
where {li)}"^ * s an orthonormal basis of Hi. Let 

n-1 

V = ( M i ® l* + 1><*1 + Mi ® + 1|) , 

'4 = 

where we assume that the operators {Mi}i satisfy certain proportionality constraints. Consider 
the Hamiltonian H = Hq + XV . Up to a constant energy shift, the effective Hamiltonian 
(cf. Section \47l\ ) is 

H eff ^(-l) n - 1 \ n (r M ---M n _ 1 r + h.c.) + O(\ n+1 ) . (6) 

The proof of this result is based on Bloch's formulation [3] of degenerate perturbation 
theory and involves a certain diagrammatic formalism adapted to the problem at hand. The 
proportionality constraints in Theorem [1] hold in the case where the Mi's are unitaries and 
are thus satisfied by the operators (J3J) and © ■ Theorem [T] therefore gives a Hamiltonian 
which generates the vertex-term A(v); the same procedure allows us to obtain a plaquette- 
operator B{p). 

We also provide a generalization of Theorem [T] to several clock systems Hi^ with po- 
tentially non-commuting operators {M|*,M^}; this is needed to generate the full quantum 
double Hamiltonian Hqd. 

Outline 

The remainder of the paper is structured as follows: In Section |3l we discuss the quantum 
double Hamiltonians and their representation in terms of loop operators. We then present 
the clock-gadget and its extensions in Section [H In Section [5l we combine these results to 
obtain our simplified Hamiltonians. We conclude in Section [51 

3 Non-abelian models based on the quantum double 

The goal of this section is to derive the necessary reformulation of the operators constituting 
the quantum double Hamiltonians defined by Kitaev [ID] . We will first review these models 
(Section [3T]) ). We then give alternative expressions for the plaquette- and vertex-operators as 
products fSection l3.2[) . While the validity of this reformulation can be checked immediately, its 
meaning is not immediately obvious, nor does it provide much insight into how our techniques 
may extend to similar models. In Section[3?3]we give what can be seen as a derivation of these 
expressions (This can be skipped at first reading) . We relate plaquette- and vertex-operators 
to ribbon-operators, and show how representations of the quantum double allow us to factor 
the latter. This is the basis of our representation of the plaquette- and vertex-operators. 
Throughout, we mostly adopt Kitaev's notation jlO) . 

3.1 Definition of the Hamiltonian 

Kitaev's models [TUJ are based on the quantum double of a finite group G. Consider a lattice 
C on an orientable 2d-surface, and assign an orientation to each edge (represented by an 
arrow). We associate a quantum degree of freedom (= C' G ') with orthonormal basis states 
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Fig. 1. An edge e with endpoints V— and v+ separating two neighboring plaquettes p_|_ and p- 



{\g)e]geG to every edge e of the lattice. The Hilbert space He of the model is the tensor 
product of these spaces. We define the following operators (indexed by g, h £ G); they act on 
the degree of freedom associated with an edge e: 

L g f {e) =E,lff*X*le L a _{e) = J2 z \zg-'){z\ e 

7*(e) =|/»)</i| e T*(e) = \hr x ){hr\ . 

Suppose that e = (v—,v+), where the arrow is directed from i>_ to t> + , and that e separates 
plaquettes p_ (on the left) from plaquette p + (on the right, in the direction of the arrow). 
That is, e has the form shown in Figure [TJ We then set 

Ls(e,v_) =L a _(e), £ 9 M+) = L B + {e) 

T h (e,p_) =Tl(e), T h (e :P+ ) = T*(e) . 

The local gauge transformations {A g (v,p)} ge G and magnetic charge operators {B g (v,p)} ge G 
associated to a vertex v on the boundary of a plaquette p are 

A g {v,p) = A g {v) = Y[ e ^ {v) L3{e,v) 

where star(u) denotes the set of edges incident to v, and where So, ... , ik-i are the edges on 
the boundary of p in clockwise order starting at the vertex v. The Hamiltonian is expressed 
in terms of the vertex- and plaquette-operators 

1 1 sec 
B(p) = B 1 (v,p) , 

where 1 £ G denotes the identity element; these are projections for every vertex v and 
plaquette p. Note that B{p) does not depend on which vertex v is chosen as the starting 
point. The quantum double Hamiltonian on T-Lc is defined as 

v p 

3.2 Rewriting plaquette- and vertex-operators as products 

In this section, we present the main property of plaquette- and vertex-operators we need. We 
show that these operators can be written as ordered products of 2-local operators acting on 
each edge of the plaquette/ vertex and an auxiliary system, followed by a projection applied 
to the auxiliary system. 
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More precisely, let Hrv = C' G be an auxiliary system placed at the center of the pla- 
quette p. The operator B(p) is a product Mq ■ ■ ■ Mf l _ l of operators Mf each acting on one 
edge e, G dp of the plaquette and the auxiliary system Hrv; the auxiliary system has to be 
projected onto a certain state \^ p )rp. A similar statement holds for the vertex operator A(v). 
The support of the operators {Mf}i and {Mf}j is diagrammatically shown in (fTJ. 

Proposition 1. Let {eo, . . . , ek-i} = star{v) be the edges incident to the vertex v, enumerated 
in clockwise order starting from eo. Similarly, let {eo, . . . , £j-i} — dp be the boundary edges 
of the plaquette p in clockwise order around p. Let Hrv = Hrp — C' G ' be auxiliary systems 
with orthonormal bases {|g)} g6 G, and let 

Tl = \g){g\ R , L a + =Y,\9z){z\R, 

z 

be projection onto \g) and left-multiplication by g on these systems. Furthermore, let 

V l G l gee 

where 1 £ G is the identity element. Define the unitary operators 

Mf = M(e h v) =Y J L 9 {e i ,v) c <8 (T 9 + ) R , Mf = Mfe.p) =Y J T a {e i ,p)c ® (^)n* • 

9 S 

on ®Hr, which act non-trivially on one edge and the respective auxiliary system. Then 

a(v) = . . . mu\^) = (^im^y ■ ■ ■ (M °)t|*") 

S(p) = (*f|M p • ■ • MjL^) = (* p |(Mj_ 1 ) t • ■ ■ (M p )t|*p) . 



Proof: This immediately follows from the definitions. □ 

To generate the complete quantum double Hamiltonian Hqd on Hc> we will later consider 
all plaqucttcs and vertices simultaneously. That is, we introduce auxiliary systems Hr« = 
Hrp = C' G ' associated with each vertex v and plaquette p, and consider the operators {Mf} Vf i 
and {Mj} p j defined by Proposition [T] These act on the joint space He <8> (&) v Hrv) <S> 
i® p H.Rv). For later reference, we show that they commute unless they act on the same edge 
of C, and are associated with a vertex and a plaquette, respectively. 

Lemma 1. Let {Mf,Mf} be defined as above. Then 

1. [Mf,M'f} = [Mf, (Mj')t] = for all v, v' and i, j. Similarly [Mf , Mf'} = [Mf, (Aff')+] = 
for all p,p' and i,j. 

2. [Mf ,Mf] — [Mf , (M?)^] = whenever these operators do not act on the same edge. 

3. If Mf = M(e,v±) and Mf = M{e,p±) act on the same edge e (cf. Figure^, then 
[Mf,Mf]^0. 
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Fig. 2. The two types of triangles corresponding to the shortest possible ribbons: the first one is 
described by an edge e and a plaquette p; it will appear in the decomposition of B(p). The second 
one is specified by a pair (e, v), and will be used to decompose the vertex operator A(v). 



Proof: To prove the first claim, we can clearly restrict our attention to neighboring pairs 
of vertices (or plaquettes p_ and p + ) that are the endpoints of an edge e = (v—,v+) 

(cf. Figure [TJ. We can further assume that = M(e, u_) and Mj + = M(e,v + ), i.e., 

that M i ~ and Mj + both act on this edge (e = e\~ = ej + ). All other commutators vanish 
because the operators act on distinct systems. Observe that the action of M i ~ on the degree 
of freedom associated with e is by right- multiplication by a group element, whereas Mj + acts 
by left-multiplication. The commutativity of these two operators therefore follows from the 
fact that left- and right-multiplication commute, i.e., [L a _{e), L+(e)] — for all g,h 6 G. This 
gives the first part of claim [TJ The second part of claim [TJ is derived in a similar fashion, using 
[T 9 _ (e),T£(e)] = for all g,heG. 

Claim [2] immediately follows from the definitions since the corresponding operators act on 
distinct systems. Finally, claim [3] can be verified by a straightforward computation using the 
definition of these operators. □ 

In the next section, we connect the expressions given in Proposition [TJ to ribbon-operators 
which create localized excitations. While this is not essential for our main result, we hope 
that this additional detail might provide some guide as to how to generalize to other models. 

3.3 Ribbon- operators, the quantum double and products 

We first recall the definition and some basic properties of ribbon-operators following [10] . 
We will subsequently use a representation of the quantum double to rewrite plaquette- and 
vertex-operators in the desired form. 

Quasiparticles in the quantum double model are associated with pairs (v,p), where v is a 
vertex on the boundary of plaquette p (these are called sites). A ribbon t connects two sites 
(v ,p) and (v' ,p'), and is defined by a path on C connecting v and v' , and an associated path 
in the dual lattice C connecting the center of p with that of p' . For every ribbon t, there 
are \G\ 2 ribbon-operators {F^ h,9 \t)} ^ . g )eG 2 which create particle pairs at the ends of the 
ribbon. These operators generate an algebra by 

F m (t)F»(t) = K n F k (t) where A ^Uh ugi ) = Shohuh 6 g0ig 6 gug . (8) 

k 

Ribbon-operators corresponding to a "long" ribbon t — toti can be written as a linear com- 
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Fig. 3. In Lemma[2] we show how the vertex-operator can be written as a sum of ribbon-operators 
corresponding to a ribbon going around vertex v. An analogous statement holds for the plaquette- 
operator B(p). This figure schematically shows these ribbons. 

bination of ribbon-operators corresponding to the ribbons to and t\ by 

F k (t t 1 ) = J2^nF m (to)F-(t 1 ) , (9) 

m.n 

where 

^\holo),(hi,gi) = 5 9,9ogi S ho < h 5 h 1 ,g- 1 hg ( 10 ) 

(Mathematically, © represents comultiplication in the Hopf algebra generated by the ribbon- 
operators, cf. [10].) The ribbon-operators corresponding to the shortest possible ribbons are 
those associated with a single triangle. These come in two types, see Figure HJ The first 
type of triangle has two sides pointing to the center of a plaquette p and an edge e on the 
boundary of p as its third side (we simply write (e,p) instead of t); the second type has two 
sides pointing to a vertex v and one side which is an edge e' of the dual lattice CJ . This type 
is specified by the edge e crossed by e' and its endpoint v. The associated operators are 

F (h ' 9 \e,p) =T 9 ~\e, P ) F (h ' 9 \e,v) = 6 g<1 L h (e,v) . (11) 

Observe that ([9]) and ([TT]) completely determine the ribbon-operators. 

We arc interested in the operators associated to closed ribbons. More precisely, we consider 
closed ribbons around vertices of the primal lattice, and similarly for the dual lattice (see 
Figure El). Let us define £[u,eo] to be the closed ribbon going around vertex v starting from 
the edge eo € star(w) in a clockwise fashion. Similarly, let t[p, eo] be the closed ribbon going 
around the center of p starting from the edge eo in a clockwise fashion; notice that these 
are products of only one type of "simple" operators. In particular, if v has incident vertices 
eo, . . . , e.j-i € star(w) (going clockwise starting from eo), and plaquette p is surrounded by 
eo, . . . ,e/c_i as explained in Section [31 these ribbons can be decomposed into a sequence of 
triangles as 

t[v,e ] = (e ,v)(e 1 ,v) ■ ■ ■ (ej_i,w) t\p,e ] = (e ,p)(ei,p) ■ ■ ■ (e k -i,p) ■ (12) 

When the choice of the starting edge eo (or eo) is irrelevant, we will drop this in our notation 
and simply write t v — t[v, eo] and t p — t\p, eo]. A straightforward computation then shows 
the following: 
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Lemma 2. The local gauge transformations {A g (v,p)} ge G an d magnetic charge operators 
{B g (v,p)}g£G ar £ ribbon- operators corresponding to closed loops, that is, 

A g (v,p) = A g (v) = F^(t v ) 
B g (v,p)=F^ 3 ' X \t\p,e }) for all h 

In particular, the vertex- and plaquette- operators can be represented 

flM = |ij e k F k (t p ) = (t p ) for all h G G 



. G G . 

as 



k 

where 



e^.ff) = s gA for all (h,g) G G 2 . (13) 



Proof: Consider a ribbon t = tot±. By ([9j, we have 

F {h ^{t t x ) = F^^^F^o'hgo^^y 

S03l=3 

Consider the case of the local gauge transformations. We show inductively that for every 
i = - 1 

F^(t ---t i )=5 g , 1 L h (e ,v)---L h (e i ,v) , (14) 
where U = (ei,v). For the base case, we get from (fTTj) 

F^(toh)= ]T 5 go . 1 L h (e Q: v)S gul LSo 1 ^o {ei7v) 

30 31=3 

= S gtl L h (e ,v)L h (e u v) . 
Assume that fj 14|) holds for i. Then 

F^\t ---UU+i)= ^• So) (*o---ti)^ (ffo " 1 ' l90 ' Sl3 (*i+i) 

3031=3 

= E **,i(n ik (er.»))«»,i £((b " I * S0,,l) (ei+i.'') ) 

3031=3 \r=l / 

where we used the induction hypothesis and (fTT|) . This immediately implies (fT4|) for all i. 
This concludes the proof of the statement about the local gauge transformations and the 
vertex-operators A(v). 

The other claims can be proved similarly, with (fT4|) replaced by 

F (h ^(t ---t l )= T<*{eo,p)-->T"(e i ,p) for all h G G , 

Si---3o=3 
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where f, = (e i: p). □ 

Our next goal is to write the plaquette- and vertex-operators as products of operators. 
We use Drinfeld's quantum double [I] T> = T>{G) of the group G. This is a Hopf algebra with 
generators {£'(/ l ,g)}(/ ! ., g )gG 2 satisfying the multiplication rule 

D m D a = ^mnAc • (15) 

k 

Note that T> is the dual of the algebra generated by the ribbon-operators. Note also that the 
local gauge transformations and magnetic charge operators ([7]) form a representation of T> by 
D(h,g) = Bh(v,p)A g (v,p). The algebra V gives a way of writing ((9]) as a product of operators. 
That is, suppose we have a representation of the algebra V. We then have the identity 

(y^F m {t a )®D m \ y^2/F n {t\) <g> Dnj =J2 Fm ( t o)F n (h)(E>D m D n 

\ m / \ n / m,n 

m Yl f4n^ m (MF n (<i)® Ac 

m,n,k 

i^F k (Wi)®D k (16) 

~ k 

In particular, by induction, we obtain the formula 

Mo--- M n _i = ^2 F k (t) ® -Dk , (17) 

k 

for a ribbon t = to • • • in— lj where 

Mi = M[U) = Yj ® D " ■ ( 18 ) 

Identity (JT7J) is our starting point for writing ribbon-operators as products. Because of 
Lemma [2j we are especially interested in certain sums of ribbon-operators corresponding to 
closed ribbons. Such a sum can be obtained from (|Tg|) by choosing a particular representation 
of T>, and projecting the second system onto a certain vector. 

Concretely, we represent the algebra I? on a Hilbert space %r = (C' G ' )® 2 with orthonormal 
basis {|k) = \k )\ki)} keG 2. The action of the generating operators on these vectors is given 
by 

Dj|k)=^n^|m). (19) 

m 

The fact that this is a representation of the algebra V follows from the identity (cf. [TDJ (31)]) 

E em ^j=E^m em =# ( 20 ) 
m m 

We define the normalized vector 



R. Konig 13 



Focusing on the case where the ribbon goes around a vertex of the primal/dual lattice, we 
get the following statement. 

Lemma 3. Let t = to • • • be of the form t v or t p (|12j) . and let Mi = M(ti) be defined 
by ([15]), Let |*) be given by (JUJ). Then 



(*|M • • • M n _x|*) = <*|M,U • • • = 



if t = t v surrounds vertex v 
B{p) if t = t p surrounds plaquette p . 



Proof: Note that the coefficients e k are real. Therefore 

mi*> = ^ el " ek 4j (22) 



|G| . 

1 1 m.k 



With (gnj), this reduces to 



<*w)= |GIV -. IGI . 



Applying this identity to (|17p using linearity gives 



if t surrounds vertex v 
B(p) if t surrounds center of plaquette p 



because of Lemma [2j Since the rhs. of this equation is Hermitian, the claim follows by taking 
the adjoint. □ 

To relate Lemma [3] to the analogous statement of Proposition [IJ we derive more explicit 
expressions for the operators Mi — M{ti) defined by (fT8|) . For this purpose, consider the 
action of the operators Dr^g)- We have 

D {hl , gi) \(h,g))= £ ^l )l(/ll , fl0 l(V<7o)> 

{ho, go) 

= J2 5 90 , gg ^^o,h5 hug -i hgo \{h Q ,g )) 

(ho, go) 

= 5 h u { ag -^h( gg -^99l X )) ■ 

From this, it is clear that 

D (h 1>gi ) = C- 1 (\h 1 )(h 1 \®l ciGl )C^(l cial ®L 9 2) , 

where C<_ is the unitary operator defined by 

C<-\(h,g)) = \{g~ x hg,g)) for all (h,g) , (24) 

and where L 9 _ denotes right-multiplication by g^ 1 . For example, if ti = (ej, v) is a triangle in 
a ribbon going around the vertex v, we get 

M(e t ,v) =J2L h (e,v)<S>C~ 1 {\h)(h\®l aG i)C^ (25) 
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by (jlip . Similarly, if ij = (e,,p) is associated with a triangle which is part of a ribbon going 
around the center of plaquette p, we have 

M{e l , P ) = Y J T 9 ~\e, P )®l c]G] ® L g _ . (26) 
g 

Observe that the operators ([25]) and ((26)) only act non-trivially on one factor of the auxil- 
iary system H.r = (C' G ')® 2 . With (j2~Tj) . (|25|) and (l26l) . the expressions given in Lemma [3] 
immediately reduce to the simpler expressions given in Proposition [TJ 

4 Perturbation theory for effective Hamiltonians 

The outline of this section is as follows: In Section 14. 1[ we give the relevant expressions 
of Bloch's perturbation expansion of effective Hamiltonians, adapted to the special case of 
interest for simplicity (see [7] for more details - we closely follow this reference) . In Section l4T2l 
we introduce our main tool, the clock-gadget. We then present an important generalization 
in Section l4~3l 

4-1 Block's perturbation series 

Let Hq be a Hamiltonian, and let Pq = J2i=i Iv^) (Vi^ I De the projection onto the d-fold 
degenerate ground space. Let {-Ey }j>o be the (ordered) eigenvalues of Hq. We assume 
that the ground state energy of Hq is Eq ^ = 0. (If this is not the case, we simply consider 
the shifted Hamiltonian Hq = Hq — Eq ^ ■ 1 since this will not affect our results.) Consider 
the perturbed Hamiltonian H = Hq + XV; let {\fi)}f =1 be mutually orthogonal eigenstates 
of H corresponding to the d lowest energies, and let {Ei]i be their energies. The effective 
Hamiltonian H c tf for the d lowest energy levels is defined as 

d 

i=l 

The perturbation XV generally leads to a shift of the energy of the entire space; we are usually 
not interested in this aspect, but only in how the Hamiltonian acts on the span of the states 
\(fii). Therefore, it is useful to introduce the shifted effective Hamiltonian 

H eS (5) =H eG -S-U, 

where n = Y2i=i ^ s the projection onto the support of H c g, and 5 parametrizes the 

magnitude of the energy shift. 

Bloch's perturbation expansion gives an expression of the effective Hamiltonian as 

iJ cff = UAU f , (27) 

with associated power series expansions of A and U in A: 

oo oo 

A =^2 A m „4 (rn) U = ^ X m U {m) . (28) 

m— rn— 
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The operator U satisfies 



(29) 



The m-th terms in the power series expansions are 

A (0) = o u (o) = Pq ( 30 ) 

and for m > 

A {m) = P VS ll VS e2 ■■■VS^- 1 VP (31) 

(d,..,<m-l)ePm-l 

w (m) = ^ S^FS^F—FS^FPo , (32) 

(«i ; ...,« m )6?,„ 

where the "reduced resolvent" is 

5 , = (E,, (-^ 0) )-^ i">o (33) 

|-P if € = 

and P m is the set of m-tuples of nonnegative integers defined as 

{m p ~| 

{li,...,£ m )\J2 £i = m >J2 £i ^ p forallp=l,...,m-ll . (34) 
i=l i=l ) 

The convergence of the series (|2"51) can be analyzed by bounding the norm of U using the 
triangle inequality (The convergence of the series for A in (|28[) follows from that of the series 
for U because of the identity A = XPqVU.), i.e., 

oo 

\\U\\ < ^2 X m \\U {m) \\ . (35) 

m=0 

This sum can further be bounded by inserting (|32[) 

||W (m) || < \V m \ ■ max \\S tl VS ta V'--VS tm VM 

<4 m max \\S ll VS e W ---VS^VPoW . (36) 
(,ti,...,i m )ep m 

Applying the submultiplicativity property \\AB\\ < \\A\\ ■ \\B\\ and the bound ||5|| < (e { ^ 
immediately gives 

(\ m 
I ■ (37) 
Ei J 

Reinserting (|37|) and (|36| into (|35[) leads to the well-known sufficient condition [7] 

E (0) 

A< 1W < 38 > 

for convergence. In our case, however, we will obtain a more refined condition by directly 
bounding the expression in (|36[) . 
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4-2 The clock-gadget: Perturbative gadgets for ordered products 

In this section, we show how to obtain an ordered product Mq ■ ■ ■ M n -\ of operators as the 
low-energy effective Hamiltonian of a 2-local Hamiltonian. Our main result is the following 
statement. 

Theorem 1. Let Hq be a Hamiltonian on a bipartite system Hs <S> Hi of the form Hq = 
—Po, where Pq — Tq <£> |0)(0|j is the projection onto the ground space. Here {l*)}"^ 1 is an 
orthonormal basis of Hi = C™. We identify \n) = |0). Let 

i=0 

where we assume that the operators {Mi}i satisfy the following proportionality constraints (we 
write A oc B if there exists a scalar c such that A = cB ): 

(i) MiM i+l Ml +l Ml oc (M,M/) 2 

(ii) (MiMjf oc M.M] 
MiMj oc M^M^ 

(w) T M a4 r Q oc r 

for all i. Consider the Hamiltonian H — Hq + XV. Let H e g(A) be the shifted effective 
Hamiltonian as explained in Section \4-. 1\ There exists a function /(A) such that for 

X < lBmaxillMill (39) 

the effective Hamiltonian is 

Heff(fW) = (-l) n ^^H target ® |0)(0|/ +0(A" +1 ) , 

where the target Hamiltonian is defined by 

Htarget = IWo ' ' ' M n - X T Q + TqM^ • ■ ■ M^T . (40) 

We point out that (|i])-([ivl) are trivially satisfied if every Mj is a unitary. 

We stress that the convergence condition (f39|) does not depend on n. This is in sharp 
contrast to the simple sufficient condition (|38|) which generally requires the gap to scale 
with the system size. In particular, we can perturbatively obtain n-local interactions from 
2-local ones without such unfavorable scaling. However, this comes at the cost of using an 
n-dimensional clock system. 

The remainder of this section is devoted to the derivation of this theorem. Since we can 
neglect constant energy shifts, we will consider the shifted Hamiltonian 



#o — Is/ — Po 
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instead of — Pq (as in the theorem) in order to be able to apply the formulas of Section |4~T1 
We first reorganize the expression pip . The projection onto the complement of the ground 
space (the (E — l)-eigenspace) is 

p 1 = p ± = i S / - r ® |o)(G|i - (i s - r ) ® |o)(o|/ + i s ® (1/ - |o)(o|/) (41) 

The reduced resolvent-operators S are given by (cf. ([55])) 

s i = {(-l)'Pi i">0 

\ -P if £ = . 

Up to a sign, the operators PqV S ll V S i2 ■ ■ ■ VS tm ~ 1 VPo therefore only depend on the set of 
indices i for which £j is non-zero. In particular, we can rewrite (|31|) as 

A im) = g m (e)P Q VP £l VP £2 ■ ■ ■ VP Sm _,VP , 

e=(ei,..., £m _i)e{0,l} m - 1 

where 

9m(e) = (-I)™" 1 J2 (-l)V<toW . (42) 

(e u ...,i m -i)ev m -i 

{i|^^0}={i|ej/0} 

Let us focus on a term of the form PoVP ei VP e2 ■ ■ - V P tm _ 1 V Pq. We define the opera- 
tors W} = M\ <E) \i + We can then write the perturbation as 

T 

where we index the family of operators Ui{Wi, W}} = {Z T } T by t. The expression of interest 
takes the form 

P VP £1 VP £2 ■ ■ ■ VP^VPo = ]T P Z Tl P ei Z T2 P S2 ■ ■ ■ P em _ x Z rm P Q . 

A similar computation can be performed for the operators Li^ m K Let us summarize what we 
obtained so far: 

Lemma 4. The operator A [m) (cf. ([31]) ) is a linear combination 

where the sum is over alls — (e\, . . . , e m — l) € {0, l} m_1 , Y = (Yi, 
the function g m is given by (|42| and where 

0(e, Y) = P Y x P ei Y 2 ■ ■ • P £ro _ 1 Y m P 

Similarly, we have 

wW=^ fa+1 (£)r( £) Y) , 



(43) 



(44) 



(45) 
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with the sum over all e = (ei, . . . ,£ m ) G {0, l} m , Y = (Yy, . . . ,Y m ) G fu* {w,-, W/j") and 

r(e, Y) = P ei YyP e2 Y-2 ■ ■ ■ Y m _yP £m Y m Po . 

Note that the operators 6(e, Y) constituting A^ m ' contain m factors Yf, we will refer to 
this as operator of order m. Our main technical result is a characterization of the operators 
0(e, Y) of order m < n. Two operators of order n play a special role; these are 

ei = e((x_^i),(w ,...,w n ^)) e: = e((L__i),(^_ 1 ,...,^ t )) . (46) 

n—l n—l 

Lemma 5. Let 0^,©£ be defined by (|46p. 

(a) All operators Q(s, Y) ^ {©i, ©Jj} of order m < n satisfy 0(e, Y) oc Pq. 

(b) ei = r M • • • Af n _!r ® |o)(o|/ ^ et = r Af+_ x • • • m%t ® |o)<o| z . 

For the proof of Lemma [5J observe that W\ increases the counter variable I from i to i + 1, 
whereas decreases it from i + 1 to i. In particular, all terms that do not match the zeros 
in the sequence (0, £i, . . . , e m _i, 0) vanish. For example, we have (for n = 5) 

e(;Wi)=P WiPo = Q 

©(•, Wl) = P W?P a = that is 0(-, Y) = for all Y = Y x 
6(11, Y)=0 for all Y= (Yi,Y 2) Y 3 ) 
J2 ©(l, Y) = PqVPxVPo = PoWoPiW^Po + fijWfUiWn-ifl) 

Y 

0(1111, Y) = PoVPxVPiVPyVPxVPo = PqWIPiWIPiWIPiWIPiWIPq + h.c. . 

Y 

We represent an operator 0(e, Y) of order m (cf. (dU)) by a diagram as follows. We first 
encode the sequence §£\ . . .e m _i0 by circles placed on the horizontal axis: is represented 
by a gray circle, and 1 is represented by a black circle. For example, 

1 

ei... e 5 = 01100 ^TiiXiTl 

-1 ; ; i i : i i 

"6"' T" i"" 3 i"" 5 "6 

To encode the sequence Y, we begin at the point (m, 0) on the horizontal axis and draw an 
upward-pointing arrow (vector (—1,1)) from the i-th to the (i + l)-th line for an operator Wj , 
and similarly a downward-pointing arrow starting from the [i + l)-th line to the i-th line for 
Wi (vector (—1, —1))- Interpreting the vertical coordinates modulo n, we have for example 

e(i,(Wo,w^) = p w p 1 w^p = y*S 
0(1. (IVj ,.ir„ ,)) = p wl_, /', u „ ; p = ^ 

"6" T"2 
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and for n — 5 



Br 




• ••••••••••• 



e 




(47) 



By definition of the operators 0(e, Y), only diagrams 

1. whose sequence of arrows forms a continuous path which 

2. starts at (m, 0) (modulo (0,n)), 

3. ends at (0,0) (modulo (0,n)) and 

4. goes through the points (£j,0) (modulo (0,n)) for all i with = 

correspond to a non-zero operator 0(e,Y) of order m. We call such diagrams valid. In 
other words, in every valid diagram the sequence of arrows must form a continuous path 
which passes through every gray circle. For example, is the sum of the following valid 
diagrams: 



A^ 



5(1H) 



5(101) h 




12 3 4 



12 3 4 



12 3 4 



12 3 4 



Using (|4T|) . we can simplify the expressions for operators corresponding to such diagrams. 
That is, we use the fact that Pi acts trivially on the subspace Hs ® span{|i)/}^0i and acts 
as (1 — To) <£) |0)(0|j otherwise. This gives for example 

PoWW + P = roAfoAfJro ® |0)(0|/ 
PqWIPxWIPiWIPiWIPiWIPq = T mIm^M^mIm^To ® |0)(0|/ . 

Essentially substituting Wi by M,, Pi by the identity (or 1 — To) and Pq by Tq, we can use 
the same pictorial representation for the resulting product of operators {Mi}i. Special care 
needs to be taken in cases where the arrows touch the horizontal axis: If the corresponding 
circle is black, the operator Ti = 1 — To needs to be inserted because of (flTj) : if it is gray, we 
need to insert r . Again, examples are 



T M M^T 



ToMl^M^ToMoM^To = 

r M,[_ 1 M n _iriMoA^ro = 



0....1...3...3. 4 



12 3 4 

In particular, the assumptions of Theorem [T] take the form 
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/ 


\ 














/ 






\ 




/N/N 



(i) 



(iii) 

(iv) i^picx • oc r , 

where the absence of circles indicates that we are away from the horizontal axis. Let us verify 
that these rules also hold for the horizontal axis. The analog of 0, combined with (|n|) is 



where the (red and blue) circles at the endpoints are arbitrary. This directly follows from 
assumption (0) and flu]), i.e., the fact that MqM\M\mI oc (M M^) 2 oc M M^. 
Next we show an analog of (jn]) , which is 



oc i 



: oc 



where the (red) circle at the left endpoint can be arbitrary. The first identity is immediate 
from assumption fnj. The second identity follows from 

M mIt x M mIy = {M Ml) 2 T - M Ml(T M MlT ) 
= c x ■ M MlT Q + c 2 • M M f r 

where we used ([n]) and (|rv)) . Finally, (pfil takes the form 



with arbitrary endpoints. This is again an immediate consequence of assumption (full) . With 
these rules, we are ready to prove Lemma [5] 

Proof of Lemma [5j It is easy to check that by repeated application of rules (|il)- (jni|) and 
their extensions, every valid diagram corresponding to an operator 9(e,Y) ^ {©ti^nl 01 
order m < n can be reduced to the diagram in (|rv|) . This implies the first claim (jaj). The 
second claim (jb| immediately follows from the definitions. □ 

Theorem Q] is an immediate consequence of Lemma [5l 
Proof of Theorem [TJ Combining d43t with Lemma [5] gives 

A [m) oc P for all m < n 

A [n) = const ■ P + <7„((1, • • ■ , 1)) • (0^ + 6+) 

n-l 

We conclude that 

n 

J2 A m -4 (m) = /(A) ■ P + 9n((l, • . • , 1))A" (FqA/o • • • Af„_!F + IW*^ ■ ■ ■ M^T ) <8 
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By (|29|) . the first summand only shifts the energy of the effective Hamiltonian H c g = UALA^ . 
Moreover, since the above identity is already a n-th order approximation to A., it suffices 
to use the 0-th order approximation U ~ = Pq to compute H e s (cf. (|30|)) to order n. 
The expression for the effective Hamiltonian in the theorem follows because g n ((l, . . . ,1)) = 

n-l 



according to (|42|) . 

To analyze the convergence of the perturbation series, we use (|45[) derived in Lemma [4] 
and the triangle inequality to get 

||^||<^| 5m+1 ( £ )|-^||r(e,Y)|| 

e Y 

< 2 m ■ 2™max| 5m+1 (e)| • max ||r(e, Y)|| . 

e Y 

Here we have used the fact that the sum over Y can be restricted to all valid paths, since 
the operator T(s, Y) is zero otherwise. Using the submultiplicativity property of the operator 
norm, we have 

||r( £ ,Y)||<||Y 1 ||...||Y m || < (maxmax{||^||,||^ t ||}) m =: 7 m . 

With the bound |g m+ i(e)| < \V m \ < 4 m , we conclude that 

\\U im) \\ < (16 7 ) m . 

From (1551) and the fact that the bounded operators form a Banach-*-algebra, that is, 



m\\ 2 = \\w}\\ 2 = iiw^ii = 



MiMi®\i)(i\ 



= || M t M} || = || || 2 



we conclude that the bound given in the theorem is sufficient to guarantee convergence of the 
perturbation series. □ 

4-3 Extension to several clocks 

In this section, we extend Theorem [I] to a situation where we have L sets of operators 
(rg 1 , {Mf}™^ 1 ) indexed by a = 1, To every index a, we associate a target Hamil- 
tonian i?target given by the analog of (|40|) . that is, 

f^gct = TMq" • • • Af„%r + r(MZ_J ■ ■ ■ (M^T where T = J] T« (48) 

We will assume that the operators {T^ } a are commuting projections, that is 

(r«) 2 = (r«)t = T£ and [r^,r^]=0 for all a,/?. (49) 

Our aim is to construct a Hamiltonian whose low-energy effective Hamiltonian is equal to 
the sum of these target Hamiltonians. We will do so in a manner similar to Theorem [T] In 
particular, we will assume that for every fixed index a, the operators (Tq , {Af"}™^ 1 ) satisfy 
the conditions of Theorem [1] That is, they obey the rules 

(i) MfMf +1 (Mf +1 )t( M f )t cx (Mf (M«)t) 2 
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(ii) (Mf(Mf)t)2 (x M«(Mf)t 

(iii) MP (Mf )t ex (Mf_ 1 )t Mf _ i 

(iv) r«M«(M«)tr«(xr« 

for all a and i. It turns out that conditions Q~ (jrv]) are insufficient for our purposes; a 
complication arises because operators with different indices a 7^ (3 may not commute. In 
addition to the above conditions, we will therefore require the following commutation relations 
for all a 7^ fj and i: 

[M i a ,rg] = [(Mf)t J rg] = o (50) 

[Mf, (M^)tM^] = [(Mf)t, K)tA^] = (51) 

We stress that we do not require commutativity of the operators Mf in the form [Mf , Mj] = 
0; indeed, these operators will in general not commute in applications of interest. We will 
take care of this non-commutativity by inserting additional operators into the perturbation. 
It will be convenient to define the sets of indices 

X (a, i) = {f3\(3^a and 3j : [Mf, M/] ^ or [Mf, (Af/)t] ^ 0} . 

To define our unperturbed Hamiltonian i?o and the perturbation V, we will introduce L aux- 
iliary systems T-Lji ® ■ • • (g> "H/i. . As before, we assume that = C" with orthonormal basis 
{|*)<*}£=o • To keep the expressions short, we will omit identities when clear from the context. 

Theorem [lj. Let (l?P , {M"}™^ 1 ) wii/i a = 1, ...,L be a family of operators on a Hilbert 
space T-L with properties (|49|) (|5 1[) and (JTJ) — (jivj) . Lei 



n-l 



be a Hamiltonian onH® ® ■ ■ ■ ® %i L , where Hi a = C™ has orthonormal basis {|i)a}i=o ■ 
Let 

V = J2 (( M ") f ® l< + ® *(«, + Mf ® + l| a ® 0) , 

where x(a,i) = ®,aex(a i) I^K^U projects onto states whose register H^b is in state |0)/3 /or 
/3 € x(cm)- Consider the Hamiltonian H — Hq + AT^. Then there exists a function /(A) 
such that for sufficiently small A 



Heff(fW) = (-1)"-^" 



0(A" 



where the target Hamiltonians H" arget are given by (|48 



We give a proof of this extension in appendix |H1 It closely follows the proof of Theorem [TJ 
with an extended diagrammatic notation for several counters. 
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Note that we have not given a criterion for the convergence of the perturbation series 
in Theorem [TJ According to the naive condition (|38|) . a scaling of the coupling strength as 
A ~ i (for constant n) is sufficient. However, this condition is not entirely satisfactory, as L 
represents the system size. A priori, it is unclear whether a significantly better criterion 
(possibly involving properties of the operators can be found. This is an important 

open problem. 

The extension from a single counter system (Theorem [JJ) to several counter systems in 
Theorem [TJ is based on coupling non-commuting terms to certain projectors (i.e., the terms 
x(a, i)) which ensure that the corresponding clocks are not simultaneously "active" (i.e., not 
in the state |0)). This allows us to treat the individual counters independently. The same 
behavior could potentially also be achieved by introducing terms in the Hamiltonian which 
assign a high energy penalty to configurations that have more than one active clock. In fact, 
this may lead to more local terms. However, with this approach, the additional terms must 
have strong coupling because the suppression of undesired configurations is only based on the 
energy denominators in the reduced resolvent. Analyzing whether this alternative approach 
yields useful results falls into the same category of problems as the previously mentioned one. 
A major difference is that Theorem [TJ only involves one energy scale corresponding to the 
parameter A, whereas the method sketched here presumably requires at least two different 
energy scales in the Hamiltonian. 

5 Perturbative gadgets for quantum double models 

We are ready to apply the clock-gadget derived in Section [JJ to the quantum double models 
discussed in Section [3] We first show how a single plaquette- or vertex-operator can be 
obtained perturbatively (Section 15. 1| ; this is based on Proposition [TJ and Theorem [JJ In 
Section [531 we then use the more general Theorem [JJ to generate the full Hamiltonian. 

5.1 Generating a plaquette /vertex- operator 

For concreteness, we consider the case of the honeycomb lattice C (it turns out that this case 
is slightly more involved than the case of a square lattice because the degrees of the vertices 
in the primal and the dual lattice are different.) We denote the Hilbert space of the qudits 
on the lattice by He- 

5.1.1 Generating a plaquette-term 

Consider a plaquette p. We give a construction of a Hamiltonian Hq and a perturbation V p 
such that the effective Hamiltonian is proportional to the plaquette-operator B(p) on He- 
Far this purpose, we introduce two auxiliary systems associated to the plaquette p: a 
counter system Hip = C 6 with orthonormal basis {|i)}f =0 and a system Hrp = C' G ' with 
orthonormal basis {\g)} ge G- The unperturbed Hamiltonian on He <& Hrp ® Hip is defined as 

H$ = ®\0)(0\ Ir , (52) 

where \^ p )rp = To define the perturbation, let eo, e\ r . . ,es be the edges on the boundary 
of p in clockwise order and let 

Mf = M(e iiP ) = T 9 (e,,p) c <8 {L 9 + ) RP for i = 0, . . . , 5 
g 
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be the operators on He ® Hrp introduced in Proposition [TJ The perturbation is 

5 
i=0 

Now consider the Hamiltonian H p = Hq + XV P . We claim that the effective Hamiltonian is 
(up to a global energy shift /(A)) equal to 

H p cS (f(\)) = -2A 6 • B(p) c ® I^X^Ibp ® |0)<0| /P + 0(A 7 ) . (54) 

To verify this statement, observe that the operators {Mf}i are unitary (cf. Proposition [1]) 
and thus satisfy the conditions of Theorem[TJ Clearly, the projection Po = (lc ® \*f?) (^\rp) ® 
|0)(0|/ P onto the ground space oi Hq also has the required form. Identity ([54")) therefore follows 
from Theorem [1] and Proposition [TJ 

5.i.£ Generating a vertex-term 

Consider a vertex v. We construct a Hamiltonian Hq and a perturbation V™ which generates 
a term proportional to the vertex operator A(v). Clearly, we could use the same procedure as 
for plaquettes. Since v has three incident edges {eo, ei, e2j-, this would give the vertex operator 
A(v) in 3rd order perturbation theory. However, this is not suitable for our purposes. Because 
our ultimate goal is to generate the full quantum double Hamiltonian Hqd, we will instead 
show how to obtain both plaquette- and vertex-operators in the same order in perturbation 
theory, with identical constants. 

As before, we introduce auxiliary systems Hr* — C G and Hi* = C 6 associated with the 
vertex v. Let 

MV = M{e h v) = Y,L 9 {e h v) c ®{Tl) R v fori = 0,1, 2 

9 

and \^ v )r^ — —h — X) eG Is) ^ e as m Proposition [TJ The unperturbed Hamiltonian has the 

yK ?l 

same form as ([52")l. that is, 

= -l^'X^'l^ ®|0)(0|/» 

and the perturbation is 

5 

V v = ^({M^ ®\i + l){i\ IP +h.c)j where (55) 

i=0 




for i = 0,1,2 
otherwise . 



It is straightforward to prove that the operators {M"}i satisfy the requirements of Theorem[TJ 
which shows (by Proposition[TJ that the Hamiltonian H v = Hq+XV v gives rise to the effective 
Hamiltonian 



#eff(/(A)) = -2A 6 • A(v) c ® |*>(tf ® |0)<0| 7 , + 0(X 7 ) 



(56) 
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5.2 Generating all plaquette/ 'vertex- operators and Hqb 

In the previous section, we have shown how to obtain a single vertex or plaquette-term in 
6th-order perturbation theory. We now consider the problem of generating several terms 
simultaneously. Our strategy is to introduce auxiliary systems "H_r« , Wiv for every vertex- 
term A(v) and Hrp, Hip for every plaquette-term B(p) we would like to generate. We use 
the same Hamiltonians as before, that is, 

H'p = -|*)<*|ja> ® |0)(0| 7P flj = ® |Q)<0|/. . 

Similarly, we use the operators M", and MV, where the superscripts indicate the different 
vertices/plaquettes these operators are associated with. 

Writing a and (3 for arbitrary plaquettes/ vertices, it is easy to see that the operators 
{Mf } a .i and the ground state projections T% = (l c ® |0)(0|/<* obey the rules (0)- 

(fiv|) (due to the unitarity of the operators {Mf} a ^) and the commutation relations (|49|) - (|5T|) 
(because operators with distinct indices a ^ ft act non-trivially on distinct systems). To 
apply Theorem [TJ, we need to consider the commutation relations between different operators 
M": According to Lemma[U we only need to take care of pairs of vertex- and plaquette-terms 
when the vertex is on the boundary of the plaquette. These do not commute if they act on 
the same edge (cf. Figure [1]). In particular, we can apply Theorem [T] in two different ways, 
giving the following statements: 

5.2.1 Generating the plaquette-part of Hqd 

Consider the perturbations V p associated to plaquette p defined by (|53|) . The effective Hamil- 
tonian corresponding to H — J2 P Hq + A J2 P V p is given by 



HedfW) = -2A 6 • (E B W) ® (® 



Up0|O)(O| /p +o(a 7 ) 



for some function /(A). An analogous statement holds for the vertex-terms. 

5.2.2 Generating the full quantum double Hamiltonian Hqd 
For every vertex v, let 

2 5 
V v = i( M i) j ® I* + ® |0)(0| F _ (e? , ® |0)(0| /P+K , + h.c) + Y, (1 ® I* + + h - c -) , 

i=0 i=3 

where p_(e?) and p+(e") are the plaquettes separated by the edge ef (this is the edge that 
acts on). We also define 



VP = E (w)* ® 1* + ® io)(oi r _ (cf) ® io)(oi r+(cf) + &.c.) , 

i=Q 

where u_(ef) and u+(ef) are the endpoints of the edge e\. The support of these operators is 

m 

visualized in ([2]). where we have omitted a diagram of the form \/*\^ corresponding to 
the second sum in the definition of V v . 
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Consider the Hamiltonian H = £ p + Y, v H o + A (£» V v + £ p V P ) ■ Then 
fleff(/(A)) = -2A 6 • Hqd ® ^(g) |*)<*|«. ® |0)(0|/^ ® ^(g) |*)<*|bp ® |0)(0|/p^J + 0(A 7 ) . 

for some function /(A). 
6 Conclusions 

Our results show how to obtain Kitaev's quantum double-based lattice Hamiltonians as the 
low-energy effective description of Hamiltonians made of simpler and more local terms. This 
is achieved with limited overhead by exploiting the relation between the Hamiltonian and 
closed anyonic (Wilson-) loops. We believe that our techniques may extend to systems such 
as Levin and Wen's string-net models, where the resulting reduction in complexity may be 
more pronounced. A major open problem concerns the convergence of the perturbation series: 
The current analysis only guarantees convergence for a coupling strength that scales with the 
system size. 
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Appendix: Proof of Theorem [T| 

We follow the steps used in the proof of Theorem [T] Again, we consider the shifted Hamilto- 
nian 

i/ =E( 1 - r ®|0)<0|a) 
a 

which has vanishing ground state energy. As in the proof of Theorem[TJ we introduce operators 

(W?)t = (Mf)t®|i + l)(i| a ®x(a,i) • 

such that the perturbation takes the form V — J2 a ii^Wt) + ^7*)- For everv fixed a, we 
index the operators L>i{W", (W 7 "") 1 '} — {Z"} T by a parameter r and sometimes also write 
V = E«, T Z? where Z? e ^{Wf, (W?)*}. 

First observe that by the projections Pq = Tq ® |0) (0| Q onto the ground space of the 
cv-th term in the Hamiltonian commute with one another. The projection onto the ground 
space of Hq is given by Pq = J\ Pq- More generally, defining P° = 1 — Pq, the projection 
onto the eigenspace corresponding to (integer) energy E is given by 

P e= E I1 P ^ E = Q,1,...,L. (57) 

(e\...,<r L )e{0 1 l} i a 
E„ s a — E 

The reduced resolvent (|33|) is 

S e = ^{-E)- l P E for £>0. 

E=£0 

A typical term in (cf. ([31]) ) takes the form 

P VS ei V---VS e ^VP = h i (E 1 ,...,E m ^)P VP El V---VP Em _ 1 VP (58) 

Ei,...,E m -i 

for the function hi = /i£ lj ... i £ m _ 1 given by 

ht{E u ...,E m -i)= ( J] (-1)^,0 J f II • (59) 

We expand the operators in this sum further using ([57]), getting 
PoVP El V---VP Em _ 1 VP = PoVP(e 1 )V---VP(e„^ 1 )VP with (60) 

E = (Ej )j=l,..., m -l 

P(e J )=P(e],...,ef) :=*[[lf ? , 
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where the sum is over all matrices e with entries in {0, 1} satisfying ^ Q ej — Ej for all 
j = 1, . . . , m — 1. Combining (|3Tj) . (f58j) and (|60|) . we obtain 

=535^7^ fe^ 51 £ ™-i J P o^P(ei)^ • • • FP(£m-i)Wo , (61) 

£ € \ a a J 

where the sums are over all i — (£i, . . . ,£ m -i) € P m -i and all matrices e = (e™)"=i'V.'m-i 
with entries in {0, 1}. 

Inserting the decomposition of the perturbation V as a sum of operators Z", the operators 
in (foTj) arc 

P yP(£i) • • • P(e m _i)VP = J! p oZ^P(ei) ■ ■ ■ P(e m -i)Z^Po ■ (62) 

(ai,-,Om) 
(Ti,...,r m ) 

Consider a term Po^" i 1 P(ei) • • • P(e m -i)Z"™ Pg. Observe that the j-th operator Z only acts 
non-trivially onH® Hi a 3 ■ Defining Y = (Y™)™ as 



: 



Zt? if ex. = otj 

1 otherwise 



we can write Z T . = ]J a Y" —: Yj, or 

P Z^P(e 1 )---P{e m - 1 )Z^P = P [) Y 1 P(e 1 )Y 2 ---P{e m - 1 )Y m P := Q(e, Y) . (63) 
In particular, this allows us to rewrite (|62|) as 

P yP(ei)---P(£ m _i)yP = ^e(e,Y) (64) 



where the sum is over a restricted set of matrices Y = (Yf*)°j of operators. Combining ((HI 
with gives the following generalization of Lemma 01 

Lemma [4j. XTie operator is a linear combination 

-4 (m) =E-9™( £ ) e ( £ > Y ) > 

where the sum is over all matrices e — {£"}"=i''''m-i wiiii entries in {0,1}, and all Y = 
(Yf*)°- with the property that for every j, there is exactly one ctj such that 

Yf = 1 for otj . 
The operators 0(e, Y) are defined by . The function g m is given by 



where he is defined by (159 
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As in the proof of Theorem[TJ we introduce a diagrammatic notation for the operators 6(e, Y). 
In essence, we stack the diagrams of Section 14.21 introducing in addition an appropriately 
placed horizontal arrow when Y" = 1. An example is 



0(e,Y) 



for 



12 3 4 





1 











3 


-1 




1 







12 3 4 



1 1 




1 1 


1 




1 




(65) 



The corresponding operator can be read off as 

e(e, Y) = PoW^PlP^P^W^iPlP^P^W^PlP^P^W^Po . 

Note that we are only interested in a subset of operators Q(e, Y), as specified by the condition 
on Y in Lemma 01. We will express this as a rule; we demand that 

0. On each vertical, there is exactly one non- horizontal arrow. 

From now on, we will only consider diagrams obeying this rule. In fact, we will further restrict 
the set of diagrams we study. As before, only certain diagrams correspond to non-zero opera- 
tors 0(e, Y), as can be seen by inspecting the definitions. In particular, the previously estab- 
lished rules apply to every "subdiagram" defined by (e a = (e™, . . . , e^), Y° = (Yf, . . . , Y£)). 
That is, every diagram corresponding to a non-zero operator 0(e, Y) satisfies for all a 

1. the sequence of arrows defined by Y Q defines a continuous path which 

2. starts at (m, 0) (modulo (0,n)) 

3. ends at (0,0) (modulo (0,n)) 

4. goes through the points (ef ,0) (modulo (0, n)) for all i with ef = 0. 

In addition to rules ©-(El), diagrams corresponding to non-zero operators 0(e, Y) satisfy 
further conditions as a result of the horizontal arrows and the operators x(a, i). In particular, 
if the circles on the verticals defined by the endpoints of a horizontal arrow are of a different 
color, then the corresponding operator vanishes. In other words, we have the additional rule 

5. for every horizontal arrow (corresponding to Y" — 1), the circles at the horizontal 
positions defined by its endpoints must have the same color (sf_ 1 — ef). 



For example, any diagram containing 
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gives 0(e, Y) = because this operator contains a product of the form PqY-P" = PqP"Yj = 
(here we used (|50|) ). Now consider the effect of the operators x(a,i). They imply that ap- 
plying Wj* (or its adjoint) from the left to an operator can only lead to a non-zero result if 
the operator is not killed by the projections |0)(0| l a for all j3 € x(a, j). This translates into 
the rule 

6. For every diagonal arrow in a subdiagram a (i.e., Yf" ^ 1), the following holds: Every 
horizontal arrow (i.e., Yf) corresponding to a subdiagram j3 € x(a,i) in the same 
vertical lies on the horizontal axis. 

For example, (|65|) satisfies this rule if and only if 1 ^ 1)- 

We will call a diagram obeying rules ((O])-© valid and restrict our attention to such 
diagrams. Again, we can write the operators 9(e, Y) associated with valid diagrams in terms 
of the operators Mi, Ti. This is done by substituting W by M and every circle that is touched 
by Tq or T\ depending on its color. For our example, we get (assuming that the diagram is 
valid) 

e(e,Y) = rM 1 (r2 r 3)(M3_ 1 )t r gM I ti(rgri])(M 1 )tr® |o)(o|® 3 . 

We now prove a few substitution rules which we will use to simplify diagrams. These 
affect two systems (a) and (f3) (and leave the others invariant). To state the rules, we only 
depict the relevant parts of the subdiagrams (e a ,Y a ) and (e^,Y^). The colored circles in 
the following diagrams may be arbitrary. 

Lemma 6. In every valid diagram, we have the following substitution rules, 
(a) 




where we assume that the endpoint of the upper sequence of arrows is not on the horizontal 
axis (apart from that, the vertical/horizontal positions of the arrows may be arbitrary). 

(b) Arrows can be commuted past triangles situated on the horizontal axis, that is, 




where we assume that the endpoint of the upper sequence of arrows is on the horizontal 
axis. 
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Proof: We show the first rule in (jaj). Let a ((3) be the index corresponding to the upper 
(lower) sequence of arrows. Because of the commutation relations (|49| -([5T |) . we can restrict 
our attention to operators with these indices. The diagram on the left then corresponds to 
an operator of the form 

(A£A£)(Mf)t(A£A£)(Mf )t(A£A*) , 

where A^,, A£ £ {1, Tq , T"} and similarly for /?. Since we assume and that this is a 

valid diagram, we have a g x(Aj), that is, [M^,Mf] = 0. Therefore 

(M/)tA2A^(Mf )t = A£(M/)t(Mf )tA^ 
= A£(M?)t(M/)tA| 

where we used (|50|) twice to obtain the first identity and the commutativity of and M". 
Multiplying this identity from the right with A^A^ and from the left with A£A^ gives 

(A£A£)(M/)t(A£A£)(Mf)t(A£A£) = (A£A£)A£(Mf )t(M/)tA£(A£A£) 

= (A£A£)A£ (Mf)+(M/)tAS(ASA^) 

= (A£A^)(Mr)t(A^)(M/)t(A«A^) 

where we used the fact that the operators A: ' are commuting projections and (|50| . This is 
the first statement in (jaj) ; the other claims can be derived in a similar manner. 

Consider the first identity in part (|bj of the lemma. Note that the diagram on the left 
corresponds to an operator of the form 

(A2A^)M «A^(M «)t(A^)(M/)t(A^) 

where A^, £ {r^T"} and A^,A^ £ {l,rQ,rf}. Again using commutativity (in particular 
assumption (1511) ) and the fact that the operators A are projections, this can be reorganized 
into 

(A2A^)M Q (M Q )t(M/)t(A^) = (A2Ai)(M/)tM «(M «)t(A^) 

= (A2Ai)(M/)t(A2A^)M^(A^)t(ASA^) 

This implies the first substitution rule in (jb| ; again, the remaining rules are shown in a similar 
manner and we omit the proof. □ 

In analogy with (|46p . we define certain operators of order n which will give us the desired 
terms in the effective Hamiltonian. These diagrams contain one of the "diagonal" diagrams 
in (j47|) as a subdiagram. They are defined as follows: For every a, we set 

= 0(e, Y) where for all j and (3 ^ a vt f* ^ °_ 



= Q(e, Y) where for all j and (3 jt a vq ^ ' J ^ i (66) 



3 3 1' 3 

Y? = lw«_J,Yf = l 



32 Simplifying quantum double Hamiltonians using perturbative gadgets 



Our aim will be to show that the operators corresponding to diagrams that are not of the 
form (|66j) are proportional to the projection Pq. We will obtain such a generalization of 
Lemma [5] by transforming the diagram in a sequence of steps into a certain form. The 
manipulations of Lemma [6] are essential for this purpose. 

More precisely, we will define three intermediate "standard" forms 51, 52 and 53 of 
diagrams. An example of the transformations we will use is 




1 2 3 4 5 




1 2 3 4 5 ! 

original 







1 


2 


3 


4 


5 


6 


7 














/ 


M 

















1 2 3 4 5 



3 4 5 

51 



52 




Lem.[9] 



1 2 3 4 5 




53 



where we indicated the lemmas explaining these transformations (we omit the circles for 
simplicity). Note that the diagram on the far right is particularly simple; the corresponding 
operator can easily shown to be proportional to Pq using a straightforward generalization of 
the contraction rules ©-(jrvj we established in Section l3~2l Thus reducing the diagram to one 
of this form is our main technical goal. 

Our first standard form 51 requires that non-trivial moves (i.e., diagonal arrows) away 
from the horizontal axis do not mix; that is, the location of the "active" arrow may only change 
from one subdiagram to another when all endpoints are on the horizontal axis. Formally, a 
diagram has form 51 if there exist indices £x,...,£ p such that (setting £$ — 0, £ p +\ = m) 

0. All chains of arrows go through (4,0), for all i. 

1. All non- horizontal arrows between £i and 4+i are located in a single subdiagram, that 
is, there is an a such that Y? € {j^W?, (W^} and Yf = 1 for alH = t t + 1, . . . , £ i+1 
and /3 a. 



For example, diagram ([55]) does not have standard form 51, but can be related to the (51- 
)standard diagram 



: o...4...a...a.. : 4 



(67) 



-1: r 
"6""l""i 3 4 



as the following lemma shows. 



Lemma 7. Consider an operator G(e, Y) ^ U Q {0"'^, S"'^} of order m < n corresponding to 
a valid diagram (e, Y). Then 0(e,Y) = 0(e',Y') 7 where 0(e', Y') is of order m with a valid 
diagram (e', Y') of form 51. 
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Sketch of Proof: This follows by iterative application of Lemma [6] (jaj) . Without loss of 
generality, assume that the first non-horizontal arrow affects e.g., Y^ = (W 1 )''. Assume 
that the horizontal axis in the first subdiagram is reached after application of the operator 
Y/ (in our example Y/ — Wq.) Clearly, all operators (arrows) between I and m affecting 
Hip , (3 ^ a can be commuted to the right of (the first upwardpointing arrow in the fist 
subdiagram) by application of Lemma [6] (jaj) . Recursive application of the same procedure to 
the resulting diagram (in particular also the part between 0, ...,£— 1) gives the claim. □ 

We will say that a diagram has standard form 52 if it consists of a sequence of "triangles" 
on each subdiagram. That is, the diagram has standard form 51 with the additional property 

2. -£ i = 2 and there is an index a such that Y£ +1 Yg +i £ {VK a (W Q ) t , (W a ) t W a }. 

Note that (1571) already has standard form 52. With Lemma [7] and the following result, we 
can reduce any diagram of interest to a diagram which has standard form 52. 

Lemma 8. Consider an operator 0(e, Y) of order m with a valid diagram (e, Y) of form 51. 
Then 0(e, Y) oc 0(e', Y'), where the latter operator is of order m' < m and where (e', Y') is 
of form 52. 

Proof: This follows by applying an appropriate generalization of the contraction rules (fijl-lfm)) 
discussed in the proof of Theorem [T] □ 

Next we introduce an additional condition: If a diagram has standard form 52 and all 
non-horizontal arrows corresponding to a subdiagram a are next to each other, we say that 
the diagram has standard from 53. Formally, this condition can be expressed as 

3. There is a sequence of distinct ax, . . . , a s and indices ki = 0, ki, ■ ■ ■ , fc s -i, k s = m such 
that Y^ e UiW.(^i ,, ) t } for a11 kj + l<k< k j+1 . 

Again, (j67|) already has standard form 53. We also have the following statement. 

Lemma 9. Consider an operator 0(e, Y) of order m corresponding to a valid diagram (e, Y) 
of form 52. Then 0(e, Y) = 0(e', Y'), where 0(e', Y') is of order m with a valid diagram 
(e',Y') of form 53. 

Proof: This follows by moving triangles corresponding to different subdiagrams past each 
other using jb| of Lemma [5J □ 

We are ready to prove the following generalization of Lemma [5J 

Lemma [5j 1. (a) Let 0(e, Y) ^ U Q {0"'^, 0"'^} be an operator of order m < n correspond- 
ing to a valid diagram, where the latter operators are defined by (j66|) . Then 0(e, Y) oc Pq. 

(b) The n-th order operators (|66p are 

e^ = rMZ---MZ_ 1 r®\o){o\® L 
0^t = r(M«_ 1 )t...( Mo Q )tr®|o)<o|» i 
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Proof: For the proof of (jaj) , we use Lemma [71 Lemma [8] and Lemma [9] This allows us 
to reduce the diagram corresponding to 0(e, Y) to standard form S3. But the operator 
associated to such a diagram is proportional to Pq, as can be seen by repeated application of 
the generalization of rules (|I])- (frvT) used in the proof of Theorem [1] 

Statement (jbj follows by inserting the relevant definitions. □ 

This is all we need to complete the proof of Theorem [TJ. 

Proof of Theorem [T]: Observe that for every configuration e with the property described 
by (l66)l . we have 

9m(e)= £ ML •■•,!)= E (-1)^^ = (-I)" 1 " 1 

fi#0 for all i 

where we inserted (|59l) and the definition (|34l) of V m -i- Lemma [U and Lemma [5[ therefore 
give 



_4. (m) ocPo for to < n 

AW = const ■ P + (-^""^target ® |0)(0|® L 



The claim of the theorem then follows as in the proof of Theorem [TJ □ 



